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ABSTRACT 

The WMAP Q, V, and W band radial profiles of temperature deviation of 
the cosmic microwave background (CMB) were constructed for a sample of 31 
randomly selected nearby clusters of galaxies in directions of Galactic latitude 
\b\ > 30°. The profiles were compared in detail with the expected CMB Sunyaev- 
Zel'dovich effect (SZE) caused by these clusters, with the hot gas properties of 
each cluster inferred observationally by applying gas temperatures as measured 
by ASCA to isothermal /3 models of the ROSAT X-ray surface brightness profiles, 
and with the WMAP point spread function fully taken into consideration. A 
reliable overall assessment can be made using the combined (co-added) datasets 
of all 31 clusters, because (a) any remaining systematic uncertainties are low, and 
(b) the data are extremely clean (i.e. free from foreground contaminants). Both 
(a) and (b) are facts which we established by examining hundreds of random 
fields. After co-adding the 31 cluster field, it appears that WMAP detected 
the SZE in all three bands. Quantitatively, however, the observed SZE only 
accounts for about 1/4 of the expected decrement. The discrepancy represents 
too much unexplained extra flux: in the W band, the detected SZE corresponds 
on average to 5.6 times less X-ray gas mass within a 10 arcmin radius than 
the mass value given by the ROSAT (3 model. We examined critically how the 
X-ray prediction of the SZE may depend on our uncertainties in the density 
and temperature of the hot intracluster plasma, and emission by cluster radio 
sources. Although our comparison between the detected and expected SZE levels 
is subject to a margin of error, the fact remains that the average observed SZE 
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depth and profile are consistent with those of the primary CMB anisotropy, i.e. in 
principle the average WMAP temperature decrement among the 31 rich clusters 
is too shallow to accomodate any extra effect like the SZE. A unique aspect of 
this SZE investigation is that because all the data being analyzed are in the 
public domain, our work is readily open to the scrutiny of others. 

1. Introduction 

One vital test of the present cosmological paradigm is the search for scattering of the CMB by 
foreground structures such as clusters of galaxies. Such observations can provide important 
information both about clusters of galaxies as well as basic cosmological parameters like Hq. 
For the CMB, scattering arises from the Compton interaction with free electrons in the hot 
(X-ray temperature) plasma of clusters of galaxies, which removes Rayleigh- Jeans black body 
flux in the direction to a cluster and leads to an apparent decrease in the CMB temperature 
- a phenomenon known as the Sunyaev-Zel'dovich effect (SZE). By now, the degree of SZE 
is highly predictable for many clusters of galaxies, because their hot intracluster medium 
(ICM) properties are well measured by X-ray satellite missions. 

In this work wc propose to perform the just such a detailed comparison. Of course, 
among the earlier papers that involved satellite data (see section 6 on interferometric tech- 
niques), the X-ray morphology and spectra of many clusters were already investigated in 
depth, while marginal detections of the SZE by WMAP were reported for individual and 
entire ensemble of clusters (Bennett et al 2003a; Myers et al 2004; Hernandez- Monteagudo, 
C, Genova-Santos, R., & Atrio-Barandela, F. 2004; Hernandez-Monteagudo, C, & Rubino- 
Martin, J. A. 2004; Afshordi, Lin, & Sanderson 2005). Nonetheless, our work in this paper 
represents the first time in which SZE in radial profiles of fluxes with errors, rather than false 
color significance maps, are shown in three separate passbands for individual clusters within 
a large (31 member) sample, leading to an in-depth evaluation of whether the WMAP SZE 
profiles arc consistent with the hot ICM profiles and temperatures as measured by X-ray 
missions. 

2. The cluster sample and predicted SZE profiles from X-ray data 

The sample employed for our purpose is the Bonamente et al (2002) catalog of 38 
nearby clusters of galaxies (hereafter simply referred to as the Bonamente sample) located 
in directions of high Galactic latitude (|&| > 30") and low column density of Galactic neutral 
hydrogen. The rationale has to do with our original intention of using the SZE as a sensitive 
probe of any extra baryonic matter that may exist in the outskirts of clusters in the form of 
a warm gas. Since this aim led eventually to a different and much more surprising finding 
we will no longer discuss it, except to mention that there is no reason why the members of 
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the Bonamente list should constitute a biased sample in any way, as far as their utility as a 
probe of the CMB distance scale is concerned. To the contrary, their high \b\ location means 
that the WMAP data for these clusters have minimal Galactic foreground contamination 
problems. 

In order to assess whether the temperature feature recorded by WMAP at the position 

of these clusters is consistent with SZE from the X-ray emitting hot ICM, it is necessary 
to first determine the radial distribution of the hot ICM density of each cluster using the 
data of X-ray observatories. In particular, we employed the ROSAT database, because the 
ROSAT XRT has a sufficiently large ficld-of-vicw to enable background determination even 
for the very extended nearby clusters in the sample, such as Coma. For clusters without 
a 'cooling flow', the standard isothermal /3-model is fitted to the X-ray surface brightness 
lx{0) in the entire region from the X-ray centroid ^ = to the background. Specifically 
three parameters are fixed by modeling the surface brightness profile: 



where is the core radius, no is the central electron number density, and (3 is the decay index. 
To predict the SZE another cluster parameters, the temperature T of the clusters, is also 
needed, and these were taken from Bonamente et al (2002). Where possible the temperatures 
were measured by ROSAT, though ASCA data were used when a cluster's ICM is too hot for 
ROSAT's passband. With the advent of Newton and Chandra, the reported temperatures 
corroborate ASCA but are generally slightly higher than the ROSAT values. Since the SZE 
depth enhances with increasing T, the Bonamente temperatures are in fact too conservative 
for estimating this depth. 

For 'coohng fiow' clusters, the central region where cooling occurs is excluded from the 
modeling. The repercussions of such a procedure arc discussed at length in section 5, where 
high resolution XMM-Newton data will be used to demonstrate that any errors in the SZE 
associated with this analysis method arc too small to affect our final conclusions. All the 
ROSAT radial profiles of the member clusters of the Bonamente sample arc shown in Figure 
1. The best-fit /5-model parameters, together with other essential details of each cluster, are 
to be found in Table 1, where a spot check against Briel & Henry (1996) reveals reasonable 
agreement with their values of Uq — 3.15 ± 0.29 x 10~^ cm~^, 9c — 5.15 ± 0.46 arcmin, and 
P = 0.93 ± 0.04 for A1795. 




(1) 



The predicted SZE decrement as a function of angle relative to the cluster center direc- 
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tion is then given by 
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where x — hv/kTcMB with v being the mean frequency of the WMAP observing filter, UTh is 
the Thomson cross section, I is the pathlength through the cluster along our off-axis sighthne, 
and the electron density is 



with r = r(/, 9) being the radius. The integration of Eq. (2) was performed analytically by 
previous authors, resulting in an expression for ATsz(6')/Tcmb which depends on no, /3, and 
the core radius Vc — L9c where L — cz/Hq (with Hq — 71 km s~^ Mpc~^) is the distance to a 
nearby cluster, see e.g. Refregier, Spergel, & Herbig (2000). We will not repeat these earher 
calculations, except to say that upon application of our Coma cluster /3-model parameters 
to the analytical formula we obtained ATsz{0 = 0)/Tcmb ~ -590 ± 170 //K at the WMAP 
frequency of z/ = 41 GHz. This compares well with the estimate by other groups, such as the 
value of ATsz(^^ = 0)/Tcmb ~ —507 ± 92.7 /xK at 32 GHz as obtained by Herbig, Lawrence, 
& Readhead (1995). 

3. The cluster WMAP temperature profiles; effect of the point spread function 

To examine whether the WMAP mission detected temperature variation in the fields of 
the Bonamente sample we extracted 'thumbnail' images (small maps) from the WMAP first 
year database centered at the X-ray centroids of the 38 clusters, i.e. the coordinates as hsted 
in Table 1. Clusters with bright radio sources along their sightlines (Virgo, A21, A1045) 
are excluded from further analysis. In addition, Fornax, A1314, and A1361 were excluded 
due to difficulties in estimating the uncertainty in the X-ray brightness distribution. Lastly, 
Hercules was not considered because the exact location of the cluster emission could not be 
identified unambiguously. For the remaining 31 sample members radial profiles of the mean 
temperature deviation over concentric annuli are computed after removing the dipole and 
quadrupole components, and plotted in Figure 2, where the error bars for each radius interval 
reflect an antenna noise which scales as the inverse-square- root of the number of observations 
of the sky area. We restrict attention to the cosmological bands of Q, V and W, since the K 
and Ka bands are dominated by Galactic foreground emissions and absorptions (Bennett et 
al 2003a,b). 

In order to compare any spatial features seen at a cluster position with the expected 
SZE behavior, it is necessary to take into account the WMAP point spread function (PSF). 
For each of the three filters we computed the average brightness profile from 15 point sources 
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and fitted it with a Gaussian function as shown in Figure 3. It was possible to perform the 
averaging of these profiles because they all compare well with each other statistically, even 
though the sources are located at very different parts of the sky, with the highest being 
at Galactic latitude |6| ^ 74°. This indicates that the performance of WMAP in spatially 
resolving CMB temperature variations is consistent across the sky. The SZE profile of each 
cluster as inferred from Eq. (2) and Table 1 is convolved with the PSF of the appropriate 
filter. A test of the correctness of our procedure was done by adopting the /?-model in Figure 
1(c) of Myers et al (2004). After convolving this model with the W-band PSF, our results 
are plotted in Figure 4. It bears close resemblance to the corresponding profile in the Myers 
paper. 

With the assurance of the aforementioned cross-checks the expected SZE decrement for 
each cluster is seen overplotted in Figure 2, where the unperturbed 'continuum' is aligned 
with the average temperature deviation in the outermost (2° - 3°) annulus. The co-added 
model and data profiles for each filter is shown in Figure 5. In all plots the error bars for 
each radius interval reflect an antenna noise which scales as the inverse-square-root of the 
number of observations of the sky area. An analysis of the (radial) bin-to-bin variation 
of the profile of many co-added random fields reveal that the degree of relative fiuctuation 
between neighboring radial bins is in accordance with the size of the error bars as shown in 
Figure 5 for each bin. This enables us to evaluate the formal statistical significance of the 
overall detection of a broad temperature decrement feature in the composite radial profiles 
as 9(7, 4.2(7, and 2.3(7 for the Q, V, and W band respectively (see Figure 5 for more details). 
Note that the true significances are hkely to be less than those given by the three numbers, 
because systematics effects in the WMAP data were not included when we calculated them. 
Such effects will be discussed at length in the next section. 

4. Systematic CMB temperature variation from WMAP random fields 

From the graphs in Figure 2, a broad CMB temperature decrement positionally coin- 
cident with the cluster and of commensurate spatial extent as the cluster size is apparent 
for some clusters, such as Coma, A1795, A1413, A2199, A2219, and A2255. On the other 
hand, there are clearly counter examples hke A85, A1367, A1689, A2029, A3301, A3558, and 
A3571. There are also clusters that fall within a 'grey' zone where no judgment can be made 
at all. Apart from some clusters having small expected SZE amidst poor signal-to-noise, 
another key problem that weakens any verdict from an individual cluster is the ambiguity in 
the CMB temperature 'continuum' appropriate to the observation. In fact, we are dealing 
with a non-trivial systematic uncertainty which can only be understood by examining many 
randomly chosen fields across the WMAP sky with |6| > 30°. For this reason we constructed 
radial profiles for 100 of such fields. 
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A plot of the r.m.s. variation (from one random field to another) of the temperature 
difference at each radial interval is shown in Figure 6, where it is evident that systematic 
effects at the level of 0.1 mK are commonplace among the smaller radius annuli. To be 
even more specific, we display in Figure 7 the radial profiles of four random fields. It can 
be seen that degree-scale modulations are frequently present, with an amplitude of ~ 0.1 
mK, often found at the ^ = (i.e. on-axis) position. The phenomenon pertains to the 
most prominent fluctuation in the WMAP data - the primary acoustic peak - which has an 
amplitude of ATcmb ~ 0.1 (or ATcmb/^cmb ~ a few x 10~^, see Bennett et al 2003b). 
Thus, when a single cluster field exhibits some broad central 'bump' it is not necessarily the 
signature of an emission source; in fact we shall demonstrate that for our cluster sample 
as a whole the discrepancy between predicted and detected SZE profiles cannot be due to 
line-of-sight emissions. Likewise, the existence of a central trough may also be caused by 
primordial acoustic oscillations rather than the SZE. The only way of knowing whether 
there is consistency between the WMAP data and cluster SZE's is to examine the co-added 
profiles of all the clusters, when the systematic effects which prevent one from determining 
the 'continuum' level can largely be suppressed. 

In order to be certain that the temperature does stabilize among co-added fields we 
generated average radial profiles of 33 random fields at a time (Figure 8). If the primary 
acoustic peaks and troughs are stacked with arbitrary central alignment in this way, one 
would expect the fluctuation amplitude to be reduced from the ~ 0.1 mK value by a factor 
of ^ -\/30, to about 0.02 mK. This is completely borne out by our averaged random profiles, 
which also indicate that towards the larger (2° - 3°) annuli the CMB temperature deviations 
from the global mean value are very uniform, hovering close to zero. By the time 100 
randomly generated profiles are averaged, the data points are smooth throughout all radii 
(Figure 8). With the help of these plots (of the merged 100 fields), wherein the systematic 
variations arc no longer a problem, we could verify if the error bar on each radial bin (antenna 
random noise) is independent of other bins. The answer is yes, because the r.m.s. scatter 
among the data points is found to be consistent with the size of the error bars. 

Returning now to the averaging of 33 random fields. Figure 8, the residual systematic 
excursions in the data are important towards understanding the radial profile of the co-added 
cluster fields. Figure 5. This is because both Figures involve the stacking of a similar number 
of WMAP fields. Since the magnitude of any remaining systematic effect is, as can be seen 
in Figure 8, far less than the discrepancy between observed and expected SZE in our cluster 
sample. Figure 5, one must conclude that the apparent incomplete SZE is not due to intrinsic 
sky variations in the CMB temperature recorded by WMAP, i.e. it is a genuine anomaly that 
deserves an explanation. 
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5. Interpretation 

How can one reconcile a cosmological CMB origin with Figure 5? It is perhaps more 
reasonable to first examine whether, despite many generations of X-ray observatories mea- 
suring the X-ray properties of clusters, we are still misled by uncertainties on the hot ICM 
parameters? Then there is the question of radio source contamination also. 

Could a steepening in the slope of the hot ICM density profile beyond the radius where 
the (5-model is well constrained by ROSAT data lead to an overestimation of the SZE? We 
pause to consider an extreme scenario under which the WMAP's spatial resolution is so poor 
that the entire SZE of even a nearby cluster simply appears as a 'point sink' at zero radius. 
If this is really the case, the WMAP radial profile of the SZE would have the shape of the 
PSF, with the flux within the entire proflle being equal to the cluster SZE integrated over 
all sightlines cutting through the cluster at various 'impact parameters'. For a /3-model, 
however, this total flux is dominated by the SZE at the outer radii where the model is no 
longer well constrained by the X-ray data. The reason is that at any off-axis angle 9 the 
line-of-sight SZE has the form: 

2 ~'~2 

(4) 

At least for our sample of clusters ROSAT data fail to guide the model only at radii 9 » 9c, 
and typically at ^ > 10 arcmin in real units. 

The total X-ray predicted SZE, integrated over all values of 9, is dominated by the 
values of A7sz(^) in the range 9 » 9c, unfortunately. This is because 

POO 

^ytotal _ / Arsz(^) 2TT9d9, (5) 

Jo 

and for 9 ^ 9c the integral ~ ^-^(/^-i), which diverges at the upper integration hmit because 
the inequality (3 < 1 applies to our clusters (see Table 1). In fact, from Eq. (3) it can 
also be seen that the total cluster mass, proportional to J ne{r)r^dr, scales with the upper 
cutoff radius in exactly the same way. Thus the outcome of our analysis is that WMAP's 
spatial resolution is very poor the X-ray model can overpredict the SZE by an arbitrarily 
large amount. This point was already raised in the recent papers of Benson et al (2004) and 
Schmidt, Allen, and Fabian (2005). Speciflcally, while Schmidt, Allen, and Fabian (2005) 
reported a steepening in the slope of the /3-model (i.e. less hot ICM) towards a cluster's 
outskirts by comparing Chandra and ROSAT data, Benson et al (2004) found that provided 
the SZE is evaluated over a radius commensurate with the resolution (or beam size) of 
the instrument any difference among the predictions derived from the various sets of model 
parameters actually remains small. 



Arsz(^) = Arsz(o) 



1 + 
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In the context of this last point of Benson et al (2004), we shall demonstrate that most 
of the results we obtained in the previous sections remain robust: the WMAP resolution is 
not as pessimistic as that depicted in the extreme scenario we just considered, because the 
detected SZE profiles are much wider than the instrument PSF. In fact, much of the SZE 
within the central 0.5° radius of 'matching' with the WMAP beam is not due to the PSF 
leaking signals from large radii inwards. Instead the opposite is true: there is a loss of inner 
signals outwards which only widens the actual gap between SZE prediction and observation. 
To test these statements, we chose a pair of /3-model parameters typical to our cluster sample, 
viz. 13 — 2/?, and 9c — 2 arcmin, and truncate the full hne-of-sight integrated SZE profile, 
Eq. (4), abruptly at ^ = 10 arcmin, resulting in the following functional dependence: 

AT^^'^'^ie) = ATsz(0)[l + 0.25(^/arcmin)2]-3 for ^ < 10 arcmin; and = otherwise. (6) 

The reason for setting the cutoff at 10 arcmin is that for most of the clusters in our sample. 
Figure 1 indicates the ROSAT data can constrain the surface brightness profile out to at 
least such a radius. If, after convolution of the truncated profile with the PSF there is a 
significant reduction of the inner SZE, this would imply a severe fiux overprediction problem 
at small radii by the WMAP PSF, which spreads signals inwards from regions beyond 10 
arcmin where the /3-model is no longer so rehable. Note that in our test we did not truncate 
the /3-model self-consistently by reducing the hne-of-sight integration to reflect the projection 
effect of the cutoff at ^ < 10 arcmin. This is a separate question to be investigated below. 

The outcome of the test is shown in Figure 9, where it can be seen that within 6 < 0.5 
arcmin the reduction in SZE by our truncation procedure is f=:i 12 % for the Q band, and 
much less for the V and W bands. This means any incorrectness in the outer (questionable) 
parts of the /3-model, when coupled with the WMAP PSF spreading effect, cannot lead to 
an overprediction of the SZE within the central 0.5 arcmin radius by a factor of four to six. 
Yet such factors are necessary to explain the average discrepancy between X-ray expectation 
and WMAP observed SZE levels for the three cosmological filter passbands, sec Figure 5. 
We must also point out that in terms of PSF effects the real problem is the opposite, viz. 
the loss of inner signals to the radii beyond, as borne out by the fact that when we convolved 
AT|^*°'^(^) with the PSF and compared the total SZE over ^ < 10 arcmin with the same 
quantity obtained before convolution, the ratio of the latter to the former ranges from 2.48 
(Q band) through 5.31 (V band) to 5.61 (W band). Since the average X-ray predicted SZE 
profile is more centrally peaked than the WMAP observed profile. Figure 5, this means the 
discrepancy between ROSAT and WMAP should have been even more pronounced if the 
original intrinsic profiles before PSF convolution were compared. 

In order to assess the hne-of-sight projection effect of the outer parts of the /3-model 
on the inner SZE prediction, we may e.g. divide the central decrement ATsz(O) into two 
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parts, respectively contributions from the hot ICM within and beyond the truncation hmit 
oi — 0rp — 10 arcmin, i.e. we write 




For our typical sample parameters of /9 = 2/3 and 6c = 2 arcmin the ratio of the 2nd term 
on the right side to the first term is 14.5 %. It is clearly not of a magnitude large enough to 
explain the discrepancy of concern, which is at the 400 - 600 % level (Figure 5). The same 
conclusion may be drawn about the other slightly off- axis sightlines of the innermost 0.5° 
radius. 

The impact of uncertainties in the X-ray /5-model on the integrity of the present work 
may quantitatively be summarized into two points. Firstly, we can calculate the significance 
of the average discrepancy between the predicted and observed SZE as depicted in Figure 5, 
taking into account: (a) the random error bars shown in Figure 5 (i.e. antenna noise for each 
radial bin, see the end of section 3); (b) the systematic error per bin which is also shown in 
Figure 5 and described in section 4; (c) the uncertainty in the /3-model parameters manifested 
as errors in the model prediction for each bin, shown in Figure 5 as the vertical interval 
between the dashed and solid lines; and (d) systematic error in the /5-modcl prediction for 
the inner radial bins due to unreliability of the model at the cluster outskirts coupled with 
the WMAP PSF smearing effect, as described in this section and quantified in Figure 9. 
When all four error components are added in quadrature, we find that for the central 1 
arcmin radius the discrepancy between prediction and observation has the significance level 
of 2.73 (T, 4.65 (T, and 4.70 cr respectively for the Q, V, and W band. This level becomes 
even higher if the next few bins are also included with our calculation. 

Onto our second summary point. As was argued earlier, since the total SZE summed 
over all sightlines out to some hmiting radius 9 ^ 9c scales with 9 as the total hot ICM 
mass (out to the same radius) does, we used Figure 5 to compute an ensemble average over 
our 31 member cluster sample the ratio of the X-ray predicted to WMAP observed total 
SZE out to ^ = 10 arcmin (i.e. 9 ^ 9c ^ 2 arcmin). This ratio, which ranges from 1.49 
(Q band) through 2.21 (V band) to 2.82 (W band), indicates by how much must the X-ray 
/?-model estimate of the hot ICM mass within the radius ^ = 10 arcmin be reduced in order 
to secure a match between ROSAT and WMAP. For the W band, which is the best (cleanest) 
cosmological filter of WMAP, the answer is close to 300 %. It is somewhat surprising that 
the /3-model can lead to such large errors in the X-ray gas mass: errors which effectively 
deemed many previous X-ray observations of the hot ICM as meaningless. Note also that, 
as already explained, because the X-ray predicted SZE profile is more centrally peaked than 
the WMAP profile, problems with the leakage of flux towards the > 10 arcmin radii by the 
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WMAP PSF render the correct percentage reduction of the hot ICM mass even higher than 
300 %. 

Could a 'cooling flow' which exists in some clusters invalidate any estimate of central 
SZE based upon isothermal (3 -model fits to the ROSAT data? While remaining in the spirit 
of the foregoing discussion, viz. on the question of how the SZE predicted profile for the 
central 0.5° radius of 'matching' with the WMAP beam may be affected by the on-axis flux 
decrement, we turn to the possible role played by the cooling of cluster cores. By using 
the more accurate measurements of Chandra, Schmidt, Allen, and Fabian (2005) found that 
a straightforward ROSAT /?-model inference of ATsz(O) could lead to an overestimate of 
the quantity by less than a factor of two. Apart from re-empliasizing that the WMAP 
observations are more discrepant from our X-ray predictions than by this factor, it must also 
be stressed that the ROSAT /3-model we employed are derived by fitting only the ROSAT 
data from regions lying beyond 'coohng flow' radii. Thus, while it is true that the central 
cooling of the hot ICM causes ATsz(O) to decrease, any accompanying central peaking of the 
gas density which has not been taken into account by the /?-model has the opposite effect. 
The net outcome depends on how the product ne{r)kT{r), i.e. the gas pressure, scales with 
radius relative to the /9-model inferred without taking the 'cooling flow' into account. After 
investigating our present sample, we found that 'cooling flows' tend to maintain or raise the 
I/O parameter rather than lower it. We illustrate this point by showing in detail the situation 
of A2029, a cluster over which the discrepancy between WMAP and ROSAT is large (see 
Figure 2). For this cluster the product ne{r)kT{r) is constrained using the high resolution 
data of XMM- Newton, and its deprojected radial profile as plotted in Figure 11 is shown to 
compare closely with the /3-model of Table 1 for A2029. 

Could line-of-sight non-CMB emissions have contaminated the WMAP passbands? In 
the restricted venue of the center of the Coma cluster Herbig, Lawrence, & Readhead (1995) 
measured a SZE of -270 //K when the prediction is -500 //K. The authors attributed the 
discrepancy to radio sources located along the central sightline but not members of the clus- 
ter. If, in general, line-of-sight sources unrelated to the clusters are bright enough to affect 
the WMAP data, the phenomenon should also be present among non-cluster directions, i.e. 
one should expect the same level of contamination to exist in the WMAP random fields. 
There is however no evidence for this, because the radial profiles of the ~ 30 accumulated 
random fields reveal an average temperature deviation of 0.005 mK for the three cosmo- 
logical bands of Q,V, and W (see Figure 8), which is on par with the average asymptotic 
deviation among the co-added cluster fields. Moreover, such amounts are far less than the 
discrepancy between the predicted and observed SZE of Figure 5. 

Could the clusters themselves be a significant source of emission in the WMAP pass- 



- 11 - 



bands? There are two possibilities. Diffuse emission from cosmic ray synchrotron radiation 
and discrete radio sources. The former is in principle a possibility (e.g. Sarazin & Lieu 1998), 
and we are currently investigating it, but see the two paragraphs after next. Concerning the 
latter, one can get a good idea of the cluster radio source occurence probability by appealing 
to the Owens Valley radio interferometry survey (Bonamente et al 2005), which finds on 
average approximately one radio source per cluster with ~ 1 mJy brightness at 30 GHz. 
Since the Owens Valley sample is more distant than ours, with a mean separation of ~ 1.5 
Gpc as opposed to ~ 0.5 Gpc for our present sample, the brightness level should scale to 
10 mJy, or 10~^^ W m'^ Hz~^. Since, to reduce the SZE by invoking emission components, 
such components must account for an average CMB temperature increase of STqmb ~ 5 x 
10""^ K distributed over the area of 0.5° angular radius. 

Is the equivalent of one unresolved 10 mJy source within the same area sufficient? The 
conversion from STcmb to a change in the observed flux involves multiplying the Rayleigh- 

Jeans sky fiux 2TrkSTcMBJ^^ / by the solid angle factor Sfl/An where Sfl = n9'^ with 9 = 
0.5°. This yields an excess flux at = 30 GHz of 10~^^ W m~^ Hz~^, ten times higher than 
the contribution from cluster radio sources. Such a conclusion is applicable to the Q band, 
which measures at a frequency only slightly higher than 30 GHz. For the W band, at a 
frequency of 94 GHz, the radio source contribution is at least two orders of magnitude short, 
because such sources typically have flux spectra ~ where a > 2. 

Since any cluster emissions which may account for the less than expected SZE detection 
by WMAP, be they diffuse or discrete in form, arc invariably non-thermal in nature, with 
a flux spectrum F^y ~ or steeper, i.e. a very different spectral shape from the Rayleigh- 
Jeans F,^ ~ u"^ dependence, another test would involve checking the WMAP band ratios of 
the SZE discrepancy. Thus a cluster source is distinguished by its larger V to W band flux 
ratios, which even in the case of a power- law as shallow as F,^ ~ is a factor of three 

more than the corresponding ratio for Rayleigh- Jeans spectra. Since the observed ratio of 
V:W for any extra radiation component that may account for the discrepancy in Figure 5 
is close to the black body limit, see Figure 10, we conclude that the resolution does not lie 
with discrete cluster radio sources. 

In addition to the above, there is one more powerful test. If clusters do exhibit a 
normal SZE which is masked by self emissions unrelated to the hot ICM, one would expect 
additional fleld-to-fleld variation of the CMB temperature at a given radial bin within our 
cluster sample, because the properties of the hot ICM and non-thermal radiation vary from 
cluster to cluster. At the very minimum, the r.m.s. fluctuation must equal that of the blank 
held and the differing degree of SZE from one hot ICM to the next, with the two effects added 
in quadrature. That even this minimum is already less than the observed r.m.s. within our 
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cluster sample is depicted in Figure 12. In fact, the observed r.m.s. is at the same level as 
that of the blank fields, excluding not only non-thermal emissions, but also the SZE itself, 
unless one contrives a scenario in which the two are correlated with each other. 

Could peculiarities in the hot ICM abundance, or the gradual decline of the hot ICM 
temperature with radius he responsible? The general reasoning can always be applied, espe- 
cially at the centers of clusters where physics are more complicated, to argue that if there are 
bright and spectrally imresolved metallic lines the density of the hot ICM electrons, which 
is ordinarily determined by assuming that over most the X-ray passband the emission is a 
continuum, may drop, resulting in a reduction of the expected SZE. In reality, however, the 
effect is very small. Thus e.g. in the case of A2029, the hot ICM abundance originally used 
by us to calculate the SZE was 0.3 solar. Analysis of XMM-Newton data revealed the same 
abundance level. Figure 11, hence no change in the density of hot ICM electrons. As to the 
large scale radial temperature dechne, numerical simulation of clusters (Romeo et al 2005) 
indicates that the hot ICM temperature can be halved by the time one reaches the virial 
radius. Observationally, however, it is clear that at least out to the radius R — 5rc beyond 
which the hot ICM contribution of the SZE is negligible (see Eq. (7)), there is no evidence 
for such a temperature drop. For instance, in the case of A2029 which has a core radius of 2 
arcmin, the hot ICM temperature as measured by XMM-Newton (Figure 11) remains stable 
out to 10 arcmin. Thus the coohng of the hot ICM towards the cluster's outskirts is not the 
reason for the large discrepancy between WMAP and X-ray fluxes. 

On the role of asphericity of either the central core or the entire cluster, and dumpiness 
of the hot ICM, these questions are harder to answer quantitatively. Given that we averaged 
over 31 clusters, however, any residual asphericity corections must be very small, unless 
one contrives a scenario under which clusters have their major axes aligned along some 
preferred direction, dumpiness is a phenomenon worthy of further research, because (to 
argue heuristically) while the X-ray observations measure (nl) the SZE probes {ug). The 
two are trivially related to each other if (ng)^ = (nl), which occurs only when the gas is 
smooth. Once clumps are present, (nl) > (ue)'^, so the X-ray prediction will overestimate 
the SZE flux. Although for the hot ICM beyond the central core it is difficult to envisage, 
from pressure balance considerations, how such high temperature plasmas could exhibit a 
large degree of clumping, the situation may be different within the core itself where the 
cooling of gas and other dynamic processes often occur. It would still be surprising if the 
factor of six discrepancy reported here for the W-band could entirely be attributed to this 
effect. Nevertheless, in the absence of further observational evidence, or a sound physical 
argument, one should not exclude the scenario of hot ICM clumping as possible cause of our 
reported discrepancy. 



-13- 



Finally, one could ask if photon populations other than the CMB may exist to interact 
with the hot ICM electrons to 'refill' the SZE flux. Whether we are deahng with 'upward' 
or 'downward' scattering, the difficulty lies with insufficient seed photon flux at frequencies 
well below and above the WMAP passbands. This is the reason for the general belief that 
cluster SZE should cause a clean removal of the CMB from its original microwave passband 
of emission. 

Conclusion 

Formally, a statistically significant detection of the SZE across our entire cluster sample 
was achieved by WMAP. The level of this detection is very weak, however. Not only is the 
field-to-field variation of the CMB temperature within the sample (due to different clusters 
exerting unequal SZE on the CMB) non-existent. Figure 12, but also the measured decre- 
ments of Figure 5 are consistent with nothing beyond the usual primary CMB anisotropy, 
viz. the systematic variations in the CMB average radial profile for a commensurate num- 
ber of co-added blank (random) fields, Figure 8. Of particular concern are the data of the 
W-band, which is the passband containing the cleanest extragalactic signals (Bennett et al 
2003b). The statistical significance of an overall W-band SZE detection is < 3 cr, when the 
expected effect is far larger. Here we also note that a similar shallow decrement was seen in 
the Myers et al. (2004) paper where the best fit model had a ATq of 0.083 mK in the W 
band. This is much smaller than the predicted average decrement from our sample of 0.46 
mK from Table 1. Thus, taken at value one may even hold the opinion that there is in fact 
no strong evidence in the WMAP database for the SZE at all, when the aggregate behavior 
of all the clusters in the sample, rather than individual cases, is considered. 

Naturally the entire premise of this paper depends on the reliability of the original 
WMAP data. If there are any data analysis issues with the WMAP processing that can 
explain the extra diffuse emission seen in our SZE clusters then our findings will be obsolete. 
However, this would have implicate severely all the WMAP analysis done to date. One 
possible resolution is to look at the SZE as probed using dedicated ground based observa- 
tories. The SZE has already been detected in a large number of high-redshift clusters using 
interferometric techniques of higher resolution than the WMAP data (e.g., Carlstrom et al. 
1996; Joy et al. 2001; Reese et al. 2002; Bonamente et al. 2005; LaRoque et al. 2005). Com- 
parison of radio interferometry and X-ray data for the same clusters show that SZE-derived 
and X-ray derived masses and gas fractions are in agreement (Grego et al. 2001; LaRoque 
et al. 2005). and allows for a determination of the cosmic distance scale (Reese et al. 2002; 
Bonamente et al. 2005). There is not necessarily a conflict between our present results and 
the previous reported SZE detections for individual clusters. As can be seen from Figure 2, 
many of the clusters in our sample do exhibit the effect at ~ the anticipated level. 
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In summary, it is through the first detailed radial profile comparison between X-ray and 
microwave observations that an apparent sample-wide discrepancy between the expected and 
measured levels of SZE from some of the best known clusters of galaxies was uncovered. The 
difficulty lies with the average behavior of our randomly selected cluster sample, which could 
still be suffering from systematic yet hitherto unknown biases. Nonetheless, the average CMB 
temperature decrement is sufficiently shallow to be interpreted simply as the usual primary 
CMB anisotropy: there is no need to invoke any SZE at this stage of the WMAP analysis. 

We are grateful to an anoymous referee for comments and critique of this work. 
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Table 1. X-ray properties of the 31 clusters employed for the present SZE investigation. The /3-model 

parameters are derived from ROSAT observations. 
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Fig. 1. — Isothermal /5-model (Eq. 1) fits to the ROSAT X-ray surface brightness data of 
the Bonamente et al sample of clusters. The ROSAT mission was chosen because of its wide 
field-of-view, which allows one to determine clearly the surrounding background level even 
for the larger clusters. For non-isothermal 'cooling flow' clusters the central region where 
the phenomenon occurs is excluded from our analysis. 
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Fig. 2. — WMAP Q, V, and W band radial profiles of the CMB temperature deviation 
(from the all sky mean value) as averaged over concentric annuli centered at the positions 
of the 31 clusters employed in our analysis. The data were plotted after subtraction of the 
contribution from the CMB dipole and quadrupole anisotropy. The solid line gives the SZE 
temperature decrement profile expected from the hot ICM with properties determined by 
X-ray observations (including the ROSAT /3-model of Figure 1) and displayed in Table 1. 
Dashed lines mark the 90 % error margin on the prediction, based upon uncertainties in the 
observed hot ICM parameters. The continuum level for the predicted SZE profile is fixed by 
aligning it with the average temperature deviation in the outermost (2° - 3°) annuli. This 
is a reliable procedure because, according to the analysis of random fields the large annuli 
averages stabilize to values near zero, see Figure 8. 
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Fig. 3. — Gaussian fits to the average profile of 15 point sources, to determine the WMAP 
point spread function for the Q, V, and W bands. The best-fit model yielded a = 0.220° for 
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Fig. 4. — Our PSF for the WMAP W band (as derived from Figure 3) convolved with the 
hot ICM model p = 0.74, 9^ = 1.5 arcmin, and amplitude AT = -0.083 K. The model 

parameters are chosen for direct comparison with Figure 1(c) of Myers et al (2004), which 
displays the temperature profile after the same model was convolved with their estimate of 
the W-band PSF. It can be seen that the two curves resemble each other very closely. This 
checks the correctness of our convolution routine. 
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Fig. 5. — The average WMAP observed and predicted radial profile for the 31 clusters of our 
sample. An individual error bar for each bin depicts the random uncertainty (i.e. WMAP 
antenna noise), while the systematic l-cr error plotted on the side of the graph depicts the 
residual large scale correlated variation in the central 1° region of the 33 co-added random 
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Fig. 6. — The r.m.s. field-to-field variation of CMB temperature at a given radial interval, 
and for the Q, V, W bands of WMAP, as obtained by comparing the radial profile 100 




Fig. 7. — The WMAP Q, V, and W band radial profile (of CMB temperature deviation) 
centered at four randomly chosen positions on the sky. 
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Fig. 8. — The WMAP Q, V, and W band radial profile (of CMB temperature deviation). 
For each band the first three graphs show the average profile of 33 fields, each centered at a 
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Fig. 9. — A comparison between the SZE predicted profile of two versions of the f3 model, 
after each profile is convolved with the WMAP PSF (for the appropriate filter). Both versions 
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Fig. 10. — Radial profile of the filter ratio of the temperature discrepancy between SZE 
model and data, for the 31 co-added cluster fields of Figure 5. First plot gives V:Q ratio, 




Fig. 11. — XMM-Newton deprojected best hot ICM parameters of A2029. The plotted 
values are averaged over all three instruments MOSl, M0S2, and PN. In the bottom graph 
the sohd line represents the /3-model as fitted to the ROSAT data outside the 'cooling flow' 
region, with the density profile then extrapolated inwards whilst assuming isothermality for 
all radii at the temperature given by Table 1 (i.e. ignoring the cooling core, which is exactly 
our procedure when we derived the predicted SZE fiux from the parameters of Table 1). 
The reason for plotting the product of deprojected density and temperature is because the 
X-ray SZE prediction depends directly on this quantity. It can be seen that our method 
yields good agreement with XMM-Newton data, even though this cluster exhibits the worst 
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Fig. 12. — Cluster ficld-to-ficld CMB temperature variation as seen in the three WMAP 
passbands at each radial bin. The solid line represents the level expected from the combined 
effect of natural blank field (i.e. primary CMB) variation and the SZE contrasts among clus- 
ters. The dashed hue represents the level expected if radio sources with average brightness 
comparable to the SZE decrement are present in each cluster with no correlation between 
such sources and the properties of the hot ICM within the same cluster. The fact that both 
solid and dashed lines are above the level of the data fluctuations is evidence for the presence 
of neither the SZE, nor radio contaminations at the level that can offset the SZE. 



